function tbrillans_f(nu,b)
    use iasing_osse
    implicit none
    ! Arguments 
    real(kind=JPR4) :: tbrillans_f
    real(kind=JPR4), intent(in) :: nu,b
 
    ! Function variables
    real(kind=JPR8) :: tbrillans_d
    real(kind=JPR8) :: nu_d,b_d
 
    real(kind=JPR8) :: nu0
    real(kind=JPR8) :: scale
 
    nu_d        = dble(nu)
    b_d         = dble(b)
 
    scale       = 100.0d+00
    nu0         = nu_d*scale
    tbrillans_d = Cst_hck*nu0/log((scale*Cst_2hc2*nu0**3/b_d)+1)
    tbrillans_f = real(tbrillans_d)
 
    return
end function tbrillans_f



